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Abstract 

Lattice data for the QCD equation of state and the magnetic susceptibility computed near the 
crossover transition at zero magnetic held are used to determine the input parameters of a five 
dimensional Einstein-Maxwell-Dilaton holographic model. Once the model parameters are fixed at 
zero magnetic held, one can use this holographic construction to study the effects of a magnetic 
held on the equilibrium and transport properties of the quark-gluon plasma. In this paper we 
use this model to study the dependence of the crossover temperature with an external magnetic 
held. Our results for the pressure of the plasma and the crossover temperature are in quantitative 
agreement with current lattice data for values of the magnetic held 0 < eB < 0.3 GeV 2 , which is 
the relevant range for ultrarelativistic heavy ion collision applications. 
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I. INTRODUCTION 

Recent relativistic heavy ion collision experiments [1-4 1 have produced a strongly coupled 
quark-gluon plasma (QGP) [5] whose physical properties are currently under intense investi¬ 
gation (see [6, 7] for recent reviews). The study of the equilibrium and transport properties 
of the QGP as functions of parameters such as the temperature T, chemical potential(s), and 
(electro)magnetic helds are of great relevance for the characterization and understanding of 
this new state of QCD matter. In particular, very strong magnetic helds up to O (0.3 GeV 2 ) 
are expected to be created in the early stages of noncentral relativistic heavy ion collisions 
[8— 13] 1 and even much larger magnetic helds of O (4 GeV 2 ) may have been produced in the 
early stages of the Universe [15, 16] (see also Fig. 10 in [17]). Moreover, magnetic helds up 

1 It is not clear at the moment if the electromagnetic fields present in the early stages of heavy ion collisions 
remain strong enough to directly affect equilibrium and transport properties of the plasma produced at 
later stages. However, see [14] for a recent study on how the QGP’s electric conductivity may actually 

delay the decay of the magnetic field in the medium. 
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to O (1 MeV 2 ) are present in the interior of very dense neutron stars known as magnetars 
[18]. Therefore, the study of the effects of strong magnetic fields on the QGP has sparked 
a large amount of interest in the community in recent years [19-46] (for extensive reviews 
and other references, see for instance, [47-50]). 

Since the properties of a strongly coupled QGP cannot be reliably studied using pertur¬ 
bative techniques one has to resort to nonperturbative approaches that are valid at strong 
coupling. Interestingly enough, contrary to what happens in the case of a nonzero baryon 
chemical potential where the sign problem of the fermion determinant prevents the appli¬ 
cation of the Monte Carlo importance sampling method in lattice simulations (for a review 
see [51]), in the case of a nonzero magnetic field (at vanishing baryon chemical potential) 
standard lattice techniques may be employed to study the equilibrium properties of QCD in 
the (T, H)-plane, see for instance, [17, 52, 53]. 

Another nonperturbative method that is suited to study strongly coupled non-Abelian 
gauge theories is the holographic anti-de Sitter/conformal held theory (AdS/CFT) corre¬ 
spondence (also known as the gauge/gravity duality) [54-56]. The correspondence has been 
employed to obtain useful insights into the properties of the strongly coupled QGP, as re¬ 
cently reviewed in [57, 58]. A very attractive feature of the gauge/gravity duality is that it 
may be easily employed to compute transport coefficients of strongly coupled non-Abelian 
gauge theory plasmas (see, for instance, [59-65]), which is a challenging task to perform on 
the lattice [66]. 

A top-down holographic dual for AT = 4 super Yang-Mills theory (SYM) in the presence 
of an external constant magnetic held was proposed in [67-69] and calculations for different 
physical observables in this scenario were carried out, for instance, in [70-73]. However, 
the QGP formed in heavy ion collisions [6, 7] probes the temperature region within which 
the QCD plasma is highly nonconformal [74] (when T ~ 150 — 300 MeV). Therefore, in 
order to make contact with realistic heavy ion collision applications, one needs to develop 
holographic models that are able to capture some of the relevant aspects of the physics of 
the strongly coupled QGP near the QCD crossover [75]. One possible way to accomplish 
this within holography is to deform the boundary quantum held theory by turning on a 
dynamical scalar held in the bulk whose boundary value sources a relevant operator in the 
gauge theory. Near the boundary the scalar held approaches zero and conformal invariance 
is recovered in the ultraviolet. In the infrared, however, the holographic dual gauge theory 
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generated by such deformation behaves very differently than a conformal plasma and may 
be tuned to display some of the properties of QCD in the strong coupling regime. 

In this work we construct a nonconformal anisotropic bottom-up holographic model that 
is suited for the study of a QCD-like plasma at nonzero magnetic held and vanishing chem¬ 
ical potential(s). Our model is built up on classical nonconformal anisotropic black brane 
solutions to the Einstein-Maxwell-Dilaton (EMD) model defined with a negative cosmolog¬ 
ical constant and in the presence of an external constant magnetic held. This constitutes 
a sequel to the studies of strongly coupled nonconformal plasmas via black brane solutions 
initiated by [59, 76] in the case of Unite temperature, zero magnetic held, and vanishing 
chemical potential", which was later extended in [62, 81] and also [63] to take into account 
the presence of a nonzero baryon chemical potential at zero magnetic held 3 . This type of 
nonconformal model has been used in the last years to investigate how different observables of 
phenomenological relevance to the QGP and the physics of heavy ion collisions vary near the 
QCD crossover transition. In fact, after the original calculations in [59, 76], which included 
the evaluation of the bulk viscosity at zero baryon chemical potential and zero magnetic held 
[59], a series of other quantities were computed within this type of holographic model such 
as the heavy quark free energy [84, 85], the energy loss of highly energetic probes [86-88], 
the Debye screening mass [77], the electric conductivity [60], a large set of hrst and second 
order viscous hydrodynamic transport coefficients [61], the spectrum of quasinormal modes 
[89] and the thermal photon production rate [90]. In the context of the holographic models 
developed in [62, 81] and [63] as extensions of the original models [59, 76], taking into ac¬ 
count the presence of a nonvanishing baryon chemical potential, we mention the calculation 
of the holographic critical point in the (T,hb )~plane and the associated critical exponents 
[81], the evaluation of the holographic equation of state, the heavy quark drag force, the 
Langevin diffusion coefficients, the jet quenching parameter, the energy loss of light quarks 
and an estimate of the equilibration time in the baryon-rich strongly coupled QGP [63], the 
evaluation of the bulk viscosity [62], as well as the baryon susceptibility, baryon conductivity, 
thermal conductivity, baryon diffusion [64], and the thermal photon and clilepton production 
rates [65] at finite baryon chemical potential and zero magnetic held. Here we add one more 

2 These nonconformal solutions can also be adapted to study the vacuum properties of the gauge theory, 

as recently discussed in [77]. This was studied in detail earlier in [78-80] in the case of similar bottom-up 

models at zero and finite temperature concerning pure glue Yang-Mills theory. 

3 See also [82] for a bottom-up holographic model at finite temperature, nonzero chemical potential, and 
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zero magnetic field in the Veneziano limit [83]. 



entry to this family of nonconformal black hole solutions by taking into account, for the first 
time, the presence of a magnetic held in the nonconformal, QCD-like gauge theory. 

Our model is a bottom-up holographic setup in which the dilaton potential and the 
Maxwell-Dilaton gauge coupling are dynamically fixed in order to describe lattice data at 
zero chemical potential(s) and vanishing magnetic held, which should be contrasted with top- 
down models coming from compactihcations of known string theory solutions. Although in 
bottom-up models the holographic dual is not precisely known, the fact that these models 
may be constructed using some phenomenological input from QCD makes it possible that 
at least part of the physics of the boundary gauge held theory resembles, even at the quan¬ 
titative level, QCD in the strong coupling limit. Thus, one may regard such constructions 
as holographic effective theories that are engineered to model some specihc aspects of QCD 
phenomenology. Once the model parameters are hxed, these theories can be used to make 
predictions about observables that are currently beyond the scope of lattice calculations, 
such as most of the second order hydrodynamic coefficients [61]. 

This paper is organized as follows. In Section II we describe in detail the construction 
of our holographic model and how the dilaton potential and the Maxwell-Dilaton gauge 
coupling can be determined by lattice data for the (2 + l)-flavor lattice QCD equation of 
state and magnetic susceptibility at zero magnetic held, respectively. With the holographic 
model parameters fully specified, we proceed in Section III to obtain the holographic equation 
of state at nonzero magnetic held and present results for the temperature and magnetic 
held dependence of the entropy density and the pressure. We fold that the deconhnement 
temperature in our holographic model decreases with an increasing magnetic held, as recently 
observed on the lattice. Moreover, our model results for the pressure and the crossover 
temperature are in quantitative agreement with current lattice data up to eB < 0.3 GeV 2 , 
which is the relevant range of magnetic helds for heavy ion collisions. We present our 
conclusions in Section IV where we also point out other applications to be pursued in the 
near future using the anisotropic nonconformal holographic model developed here. 

Throughout this paper we use natural units c = h = ks = 1 and a mostly plus metric 
signature. 
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II. THE HOLOGRAPHIC MODEL 


Assuming as usual that charm quarks are not relevant in the crossover transition, in 
QCD there are three different chemical potentials associated with three independent globally 
conserved charges. These different chemical potentials are the three lighter quark chemical 
potentials fi u , /i^, /i s or, equivalently, the baryon chemical potential /i#, the electric charge 
chemical potential /tq, and the strangeness chemical potential ^s■ For each nonzero chemical 
potential in the gauge theory there must be a nonzero temporal component of the associated 
gauge field in the bulk. It is also clear that an Abelian magnetic field B in the gauge theory 
should come from a nonzero spatial component of the gauge potential in the electric charge 
sector . 1 

In the present work we solely focus on the electric charge sector at B ^ 0 with /tq = 
Bb = Bs = 0, which may be described by the following EMD action 


S = 


167tG 


d 5 x\/^g 


5 Jm 5 


R 


\(WY - V(4>) - ^f-F 2 


flV 


+ Sghy + S { 


CT, 


( 1 ) 


where S'ghy is the Gibbons-Hawking-York action [91, 92] needed to establish a well-posed 
variational problem with Dirichlet boundary condition for the metric, and Sqt is the coun¬ 
terterm action that can be constructed using the holographic renormalization procedure 
[93-97]. These two boundary terms contribute to the total on-shell action but not to the 
equations of motion and, since we shall not need to compute the total on-shell action in the 
present work, we do not need to worry about their explicit form here. Also, as we are going 
to discuss in detail in Section IID, we shall dynamically fix the gravitational constant G 5 , 
the dilaton potential V((/)), and the Maxwell-Dilaton gauge coupling /(</>), by solving the 
equations of motion for the EMD fields with the requirement that the holographic equation 
of state and magnetic susceptibility at zero magnetic field match the corresponding lattice 
QCD results. 

In (1), the metric field in the bulk is dual to the stress-energy tensor of the boundary 
field theory while the dilaton field is introduced in order to dynamically break the conformal 
symmetry of the gauge theory in the infrared. The Abelian gauge field in the bulk is 
employed here to introduce an external magnetic field at the boundary, which we take to be 


4 See, for instance, [69] for a construction with finite B and finite electric charge density in the context of 
Einstein-Maxwell-Chern-Simons theories. 
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constant and uniform in the 5-direction and, as stated before, in the present work we set all 
the chemical potentials to zero. The constant and uniform magnetic held breaks the 50(3) 
rotational invariance of the gauge theory down to 50(2) rotations around the 5-axis implying 
that the Ansatz for the bulk metric must be anisotropic and translationally invariant. Also, 
at zero temperature this Ansatz must be invariant under boosts in the (£, z)-plane though 
this symmetry is not present at nonzero temperature. Based on these symmetry properties, 
which are phenomenologically dictated by the corresponding symmetry content present in 
current lattice QCD calculations defined on the (T, 5)-plane, we take the following black 
brane Ansatz for the bulk fields in 5 (1): 


ds 2 = e 2a(r) [-h(r)dt 2 + dz 2 ] + e 2c(r \dx 2 + dy 2 ) + 


e 2b ^ r ) dr 2 
h{r) 


(f> = 0(r), A = A^dx^ = Bxdy =>■ F = cL4 = Bdx A dy , (2) 

where the radial location of the black brane horizon, r#, is given by the largest root of 
the equation h(rn) = 0 and in our coordinates the boundary of the asymptotically AdSs 
spacetime is located at r —>■ oo. In (2) we have already fixed a convenient gauge for the 
Maxwell held, which in the present case is a prescribed non-dynamical held. Also, for 
simplicity, we shall adopt units where the asymptotic AdSs radius is equal to one. 

Using (2), the equations of motion obtained from (1) may be expressed as follows 


+ 2 a' + 2c 


h'\ „ e 2b /<9U(0) B 2 e 


2 4c 


6 ' + ftU'“TV d<p 


+ 




d(p J 


= 0, 


„ (14 , ,, 4 ti\ , 8 /2 2 /2 2 h' , 2 e 2b . . 1 , n n 

a + -b+-jja+-a + -c +--c + (<f>) - ~4> - 0 , 

„ /10, 1 h'\ , 2 /2 4,2 2 ti , le 2 * 1 ,/2 n 

c - (t“ + 6 7 + 3 h) ^ + 3 C - 3° - 3l a ~ + = °- 

h" + (2a' + 2c' - b ') ti = 0, 


(3) 

(4) 

(5) 

( 6 ) 


where the prime denotes a derivative with respect to the radial direction. Using these 
equations of motions one can also derive a useful constraint 


a 


ti 


-,2b 


a' 2 + c ' 2 - T 0' 2 + ( ^ + d J £ + 4a'c' + ^ ( V{<f>) + 


2 4e 


B 2 e 


-m = 0 . 


(7) 


5 As we shall discuss soon, B is one of the two initial conditions controlling the temperature and the external 

magnetic field at the boundary quantum field theory. The other initial condition corresponds to the value 
of the dilaton field evaluated at the black brane horizon, <j> o- The set of initial conditions (</>o,£>) is 
nontrivially related to the thermodynamical pair (T, B) in the gauge theory. In Sections IIC and IID we 
discuss how one can relate B to the external magnetic field at the boundary gauge theory, B. 
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The equation of motion for the Maxwell held is automatically satisfied by the Ansatz (2). 
Moreover, b(r) has no equation of motion and, thus, it can be freely chosen to take any value 
due to reparametrization invariance. In the next Section we specify a subsidiary condition 
for b(r) that defines a convenient gauge for the metric that will be used in the numerical 
calculations carried out in the present work. 


A. Ultraviolet expansions 


For the calculation of physical observables in the gauge theory one needs to obtain the 
near-boundary, far from the horizon expansions for the bulk fields a(r), c(r), h(r), and 
4>(r). In the present work, we use the domain-wall gauge defined by the subsidiary condition 
b(r) = 0. At the boundary the dilaton held goes to zero in such a way that V(<p(r —>• oo) —)• 
0) = —12 and /(0) is a finite positive constant 6 . Also, the metric blackening factor, h(r), 
must go to a constant at the boundary, which we denote by' h(r —> oo) = . 

Moreover, since we are interested in asymptotically AdS 5 solutions to the equations of 
motion (3), (4), (5), and ( 6 ), at the boundary one hnds a(r —» oo) = c(r —>■ oo). In the 
domain-wall gauge b(r) = 0 , the leading order near-boundary expression for a(r) (and also 
c(r)) is linear in r [62, 81] such that at lowest order in <f)(r —> oo) —y 0 we may consider the 
following leading order far from the horizon ultraviolet asymptotics 


Vtf)*-12 , /(0W(0), h(r)^h^, a(r)«4 ar + a f -r, c(r) ^ (8) 

where = c f *\, as discussed above. Indeed, by substituting ( 8 ) into the equations of 
motion and taking the asymptotic limit of large r (where the ultaviolet expansions hold), 
one concludes that 


a fa i = 


1 

vW 


( 9 ) 


In order to obtain the next to leading order term for h{r) and also the hrst terms for <f>(r) 
in the ultraviolet expansions for the bulk fields, we consider the first backreaction of the 
6 Note that in (1) the Maxwell-Dilaton gauge coupling /(</>) plays the role of an inverse effective gauge 


coupling squared and, therefore, it must correspond to a positive-definite function. 

7 This constant is equal to one in the so-called “standard coordinates” of the domain-wall gauge, which we 


shall discuss soon. Here we are considering general coordinates where this constant may be different than 


one. We shall also see later how to relate these two sets of coordinates. 





near-boundary fields expressed in (8) and (9) on the equations of motion 8 . In fact, we first 
consider the next to leading order near-boundary expansion for the dilaton potential 

9 
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V{<j>) ~ —12 + —0 2 , 77(2 = —z'A, (10) 

where A is the ultraviolet scaling dimension of the gauge invariant operator dual to the bulk 
dilaton field and we defined v = d — A, where d = 4 is the dimension of the boundary. We 
shall see in Section IID that a good description of lattice data can be achieved by taking 
A ph 3 (u ~ 1). One can now show that the far from horizon ultraviolet asymptotics for the 
bulk fields may be written as 

a(r ) ~ a(r) + • • • , 
c(r)^a(r) + (4 ar -4 ar ) + -- - , 
h(r) « h { 0 ar + hf e~ 4a(r) + • • • , 

</>(r) « (j> A e-^ r) + ^e~ A “ (r) + • • • , (11) 

where we defined a(r) = aQ ar +r/ y/h^ while • • • denotes subleading terms. We note that the 
ultraviolet asymptotics (11) are in agreement with our numerical solutions. By comparing 
these numerical solutions to (11) one can determine the ultraviolet coefficients ag ar , Cg ar , h Q ar 
and <pAi which are needed to compute the thermodynamical observables in Sections IIC and 

IID. 

B. Infrared expansions 

Now we consider the infrared, near-horizon expansions for the bulk fields a(r), c(r), h(r), 
and 4>(r). Near the horizon all the bulk fields in (2) are assumed to be smooth and we may 
consider the Taylor expansions 

OO 

X(r) = ^X„(r-r H )», (12) 

n =0 

where X = {a, c, h, </>}. 

In order to numerically solve the equations of motion (3), (4), (5), and (6) we need to 
specify the boundary conditions A(r star t) an d A'^tart), where r start is a value of the radial 
8 This procedure may be repeated to obtain all the other subleading terms in the ultraviolet expansions. 
However, we only need the first few terms in these expansions to compute the thermodynamical observ¬ 
ables. 
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coordinate that is slightly above the horizon 9 . In this paper we work with Taylor expansions 
up to second order, which are sufficient to perform the numerical integrations if r star t is 
close enough to r#. Therefore, we must determine 12 Taylor coefficients in order to specify 
X (ystart.) and X'(r start) at second order. One of these 12 coefficients, namely, 0 O , is one of the 
two initial conditions of the problem 10 . Four of these 12 coefficients, namely, do, Co, ho, and 
h± and also the radial location of the black hole horizon, r#, may be fixed by rescaling the 
bulk coordinates while taking into account also the fact that h(r) vanishes at the horizon. 
For definiteness, we adopt here numerical coordinates fixed in such a way that 

r H = 0; a 0 = c 0 = h 0 = 0, hi = 1. (13) 

Note that r# = 0 may be obtained by rescaling the radial coordinate while ho = 0 comes 
from the fact that h(r ) has a simples zero at the horizon. Also, hi = 1 may be obtained by 
rescaling t while do = 0 may be arranged by rescaling (t, z) by a common factor. Similarly, 
Co = 0 may be arranged by rescaling (x, y) by a common factor. After this, the remaining 7 
coefficients in the near-horizon Taylor expansions for the bulk fields can be fixed on-shell as 
functions of the initial conditions (0o> £>) by substituting the second order Taylor expansions 
into the equations of motion and setting to zero each power of r start in the resulting algebraic 
equations 11 . 

With A (r star t) and -X 7 (rstart) determined as discussed above, the equations of motion 
are numerically integrated from r star t near the horizon up to some numerical ultraviolet 
cutoff r max near the boundary. We used r sta rt = 10 -8 and r max = 10 to numerically solve 
the equations of motion. It is important to remark, however, that even before reaching 
Aonformai = 2 the numerical backgrounds we considered in the present work have already 
reached the ultraviolet fixed point corresponding to the AdSs geometry. This fact is used in 
Section IIC to reliably obtain the ultraviolet coefficients in (11) and it will be also employed 
in Section IID to properly compute the holographic magnetic susceptibility numerically. 

9 The horizon is a singular point of the equations of motion and, thus, we need to initialize the numerical 

integrations slightly above it. 

10 As discussed before, the other initial condition is B. 

11 In practice, we set to zero the following 7 terms: 0(rg taTt 

^start ), an d fl(r s tart) (tl) ? fl(r s tart) nt (1); 

ChAIart) and (’’start) in (3), and 0(r° tart ) in (4). 
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C. Coordinate transformations and thermodynamical observables 


Let us now introduce the so-called “standard coordinates” of the domain-wall metric 
gauge, b(f) = 0, where variables with ~ refer to quantities evaluated in these standard 
coordinates where the background reads 


ds 2 

4> 


g 2a(r) _ 

4>{r), A 


h(f)dP + dz 
= Af/ix' 1 = 


2 + e 2a ( f )(dx 2 + dy 2 ) + A -—, 

-I h(r) 

Bxdy => F = dA = Bdx A dy, 


(14) 


and the boundary is at f —>• oo while the horizon is at f — fn- The “hat” in B accounts 
for the fact that this is the magnetic field measured in units of the inverse of the AdS 
radius squared, while B shall be used to denote the boundary magnetic held measured in 
physical units, as we shall discuss in Section IID. In the standard coordinates, the ultraviolet 
asymptotics for the bulk fields are given by [62, 81] (see also [63]) 


a{r) 

c(r) 


iy> j ... 
iy> j ... 


J 


J 


h(r) ~ 1 + • • • , 
4>(r) ~ e~ vr + • • ■ . 


(15) 


The standard coordinates (in which h(r) goes to one at the boundary) are the coordinates 
where we obtain standard holographic formulas for the gauge theory’s physical observables 
such as the temperature and the entropy density. However, in order to obtain numerical so¬ 
lutions for the bulk fields one needs to give numerical values for all the infrared near-horizon 
Taylor expansion coefficients, which in turn requires rescaling these standard coordinates, 
as discussed in the previous Section. The numerical solutions are obtained in the numerical 
coordinates described by the Ansatz (2) with the ultraviolet asymptotics (11), while stan¬ 
dard holographic formulas for physical observables are obtained in the standard coordinates 
described by the background (14) with the ultraviolet asymptotics (15). One may relate 
these two sets of coordinates by equating 0(f) = 0(r), ds 2 = ds 2 and Bdx A dy = Bdx A dy 
and this leads to the following relations 12 (by comparing the near-boundary asymptotics 

12 As mentioned in [81], if (f>A < 0 one must replace 4>a >—> \4>a\ in these relations. 
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(11) and (15) for r —>• oo) 


r = 


+ a|f r - In 




t = 


X = 


g c o a o x , 


y = a o y, 


^far_^ far 

_l 0 


5 


2 A / j 


a(r) = a(r) — In 

c(f) = c(r) - (4 ar - a(? r ) - In (VY") , 
h(r) 


h{r ) = 


h% 


o 


<f>{r) = 0(r); 

e 2 (^ ar -4 ar ) 


5 = 


>r 


B. 


(16) 


The temperature of the plasma is given by the black brane horizon’s Hawking temperature 


T = 


-y'rjr' 


47T 


r=r H 


oa{f H ) 

47T 


>'(^)| = 


4<#’ 


(17) 


while the entropy density is obtained via the Bekenstein-Hawking’s relation [98, 99] 


S= H = 


5 _ A h /AG 5 = /horizon d 3 x\/g(f = f H , t fixed) = 2vr a(r ^ )+2£(r - ff) = 27 re 2 ( a o ar 4 ar ) 

4G 5 H k 2 


y 


o / 3/y 

K 2 <Pa 


(18) 


where we defined k 2 = 87rGs and used (12), (13), and (16). 

One can see from (16), (17), and (18) that the only ultraviolet coefficients in the numerical 
coordinates which we need to fix by fitting the numerical solutions with (11) are OQ ar , Cq, 
h g ar , and (j) A ■ The numerical solutions for h(r) converge quickly to their asymptotic values 
at large r and we may reliably set = h(r con formal)- With fixed in this way, we may 
fix aQ ar , 4 ar , and 4>a , respectively, by employing the fitting functions a(r) = a g ar + r/y/h g ar , 
c(r) = Cg ar + r/y/h^, and 4>{r) = (t) A e~ ua{r) in the interval r e [r con formal - 1 , r CO nformai] • We 
were able to obtain good fits for the near-boundary behavior of the numerical solutions using 
this fitting scheme. 
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Also, it is important to remark that there is an upper bound on the initial condition 
B for a given value of the initial condition for the dilaton (j) 0 . In fact, for values of B 
above this bound, all the numerical backgrounds we generated failed to be asymptotically 
AdSs. Such a bound, which we denote by B < £> max (0o), may be numerically constructed by 
interpolating a list with pairs of points {(</>q, -^max), i — 1, 2, 3, • • ■ } and the corresponding 
result is presented in Fig. 1. 



FIG. 1. The curve corresponds to the upper bound for the initial condition B as a function of 
the initial condition for the dilaton 0o, below which the solutions of the equations of motion are 
asymptotically AdSs. This curve depends on the chosen profiles for the dilaton potential V((f>) and 
gauge coupling function f(<f>) to be discussed in the next Section. 

In the next Section we explain how one can express the thermodynamical quantities B , 
T, and s in physical units 13 using the lattice data for the equation of state and the magnetic 
susceptibility at zero magnetic held. 

13 Note from (16), (17), and (18) that B, T, and s are proportional to 4>a ^, 4>a^ U and ^a^- respectively. 
Correspondingly, their counterparts in physical units (without the “hat”) are given in MeV 2 , MeV, and 
MeV 3 , respectively. This is related to the fact that the leading mode for the dilaton field, <f>A, corresponds 
to the insertion of a relevant deformation in the quantum field theory, which is responsible for generating 
an infrared scale that breaks the conformal invariance of the theory at low energies [62], 


13 





D. Fixing the Maxwell-Dilaton gauge coupling using lattice data for the magnetic 
susceptibility at zero magnetic field 


Refs. [61, 63] discussed in detail how to dynamically fix the dilaton potential, V (0), and 
the gravitational constant, k 2 , using the recent lattice data [74] for the QCD equation of 
state with (2 + 1 [-flavors. We refer the reader to those papers for the details about this 
procedure. The results are 

V(<j>) = -12 cosh(0.606 <j>) + 0.703 </> 2 - 0.1 0 4 + 0.0034 0 6 , k 2 = 8ttG 5 = 12.5 . (19) 


From the dilaton potential specified above one obtains the dilaton mass m 2 « —3, as antic¬ 
ipated in Section II A. 

We remark that, although the present EMD construction do not explicit introduce fun¬ 
damental flavors at the dual boundary quantum held theory, the dilaton potential in Eq. 
(19) was adjusted in order to quantitatively mimic the (2 + l)-havor lattice QCD equation 
of state and its crossover. This mimicking procedure was originally introduced in [76] (see 
also [100 1 for more recent discussions), where it was also discussed how different choices for 
the dilaton potential may emulate not only the QCD crossover, as done in the present work, 
but also first and second order phase transitions, which may be useful for a large variety of 
different physical systems. 

In the present paper, we employ the same procedure used in [61, 63] to express the 
holographically determined thermodynamical observables in physical units, i.e., we find the 
temperature at which our speed of sound squared, c 2 , displays a minimum (at zero magnetic 
held) and match it to the corresponding lattice QCD result [74] 


A = 


T 


lattice 


T BH 2 

min. ci 


143.8 MeV 
0.173 


831 MeV. 


( 20 ) 


In what follows, we relate any black hole thermodynamical observable, X, with its counter¬ 
part in physical units, A", with mass dimension [MeV p ], by taking X = X P X [MeV p ]. This 
prescription respects the fact that dimensionless ratios, such as s/T 3 , must give the same 
result regardless of the units. A comparison between our holographic results for the speed 
of sound squared, c 2 (T,B = 0), and the (normalized) pressure, p(T,B = 0)/T 4 (at zero 
magnetic held) and the corresponding lattice QCD results from [74] is shown in Fig. 2. One 
can see that the holographic model provides a good description of the lattice data in the 
absence of an external magnetic held. 
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T [MeV] 



FIG. 2. Holographic calculation of the speed of sound squared (? s and the (normalized) pressure 
p/T 4 . The data points correspond to lattice QCD results from [74] computed at zero magnetic 
held. 


In order to fully determine our holographic model and include the effects from a magnetic 
held we also need to dynamically fix the Maxwell-Dilaton gauge coupling /(</>). This can be 
done using the recent lattice data [52| for the magnetic susceptibility of QCD with (2 + 1)- 
flavors evaluated at zero magnetic held. In order to compute the magnetic susceptibility in 
our holographic model we follow the same general steps discussed in [101]: we substitute 
the Ansatz (2) into the action (1) and calculate the second derivative of the on-shell action 
with respect to the magnetic held, dividing the result by the entire spacetime volume of the 
boundary. In order to obtain the bare magnetic susceptibility we plug the on-shell numerical 
solutions into the expression obtained in the previous step 11 , 


Xb»re(r,B) 


ay, 


bare 


dB 2 


1 &S%g£[B] 

hbdy dB 2 


1 d*SZ? en [B] 


hbdy 

1 

“2^ 


dB 2 


^fixed 
' max 


drf(4>(r))e 2 d a ^-^d 


f r H 


on-shell 


( 21 ) 


where /b are is the bare free energy density and, formally, one should take the limit —y 

oo. However, in numerical calculations, r 1 ™) 1 must be a hxed ultraviolet cutoff for all the 
geometries in order to ensure that the ultraviolet divergence in (21) is independent of the 
temperature. Since we are interested here in calculating the magnetic susceptibility at zero 
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As mentioned in footnote 7 of [101], the Euclidean action has the opposite sign of the Lorentzian action. 
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magnetic field where a(r) = c(r), one obtains from (21) 

1 


,*;fixed 
' max 


Xbar e(T, B = 0) = 


2k 2 


drf((j)(f )) 


r H 


on-shell 


( 22 ) 


In order to regularize (22) we follow the same procedure adopted on the lattice [52] and 
subtract from (22) the vacuum contribution at zero temperature. Clearly, this removes the 
ultraviolet divergences since those are temperature independent. More precisely, we subtract 
the geometry corresponding to (T sma n, B ) « (0.005 MeV, 0), which is generated by the initial 
conditions (bo, B) = (7.8,0); this is the asymptotically AdSs geometry with the lowest 
temperature and zero magnetic field which we could reach in our numerical computations 15 . 
Therefore, we obtain the following holographic formula for the magnetic susceptibility at 
zero magnetic held (which is valid for any EMD model of the kind considered here) 

X(T, B = 0) = Xbar e{T, B = 0) - Xbare^mall, B = 0) 



affixed \ 

/ ' max \ 



/ drf(4>(f)) 

— (same) 


Jr H ) 

T,B =0 

T S ma.\l,B=0_ 


on-shell 


r var \ 

/ drf{<t>(r )) 

— (same) 


'^T’start J 

T,B =0 

T sman ,S=0_ 


I on-shell 


(23) 


where must be chosen in such a way that the upper limits of integration in the nu- 

-fixed 


r var 

max 


= VW 


- a { 0 ar + In 


A 


< r 


_ ' max 


for 


merical coordinates satisfy r con formai < r 
all the geometries considered. We found that for ~ 33 such requirement is met. We 
also checked that one can vary the value of the ultraviolet cutoff and the results for 
the holographic magnetic susceptibility do not change, which confirms the stability of our 
numerical procedure. 

We can now use many different trial profiles for /(</>) to evaluate (23) over the zero 
magnetic field background solutions, trying to holographically fit the recent lattice data 
from [52] for the magnetic susceptibility of (2 + l)-flavor QCD at zero magnetic field. We 
found that a good description of the lattice data can be obtained by fixing 


/(</>) = 1.12 sech(l. 05 0- 1.45), (24) 

15 Note that (f> o = 7.8 corresponds to the local minimum of our dilaton potential (19). For (f> o > 7.8, our 
dilaton potential becomes non-monotonic and, in practice, we took bo = 7.8 as the upper bound for the 
initial condition bo in our numerical calculations to avoid complications with extra singular points in the 
equations of motion. 
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FIG. 3. Holographic calculation of the magnetic susceptibility at zero magnetic held and comparison 
with lattice data from [52] (we consider 10.9 times the data available in table III in [52], which 
corresponds to the magnetic susceptibility in natural units - see footnote 1 in [52]). 

with the corresponding results displayed in Fig. 3. 

With the dilaton potential (19) and the Maxwell-Dilaton gauge coupling (24) dynamically 
fixed by the description of adequate lattice data at zero magnetic held, our holographic 
model is now fully determined. This setup may be employed to investigate the physics 
of the dual quantum held theory at finite temperature and nonzero magnetic held with 
vanishing chemical potential(s). 

We finish this Section by mentioning some limitations of the holographic model presented 
here: 


• The model cannot describe phenomena directly related to chiral symmetry and its 
breaking/restoration (such as T = 0 magnetic catalysis [102-104]). This could be 
studied by adding havor D-branes in the bulk (see, for instance, Ref. [82]); 


• The model cannot properly describe hadron thermodynamics (which sets in at low tem¬ 
peratures, below T ~ 150 MeV) and the effects of magnetic helds at low temperatures 
(for a study of the hadron resonance gas in a magnetic held see [105]). Moreover, 
in this holographic model asymptotic freedom is replaced by conformal invariance 
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at sufficiently high temperatures. Furthermore, for high enough magnetic fields the 
nonlinear nature of the DBI action for the D-branes should be taken into account [106]; 


• In Appendix A, we present a brief discussion on the behavior of electric field response 
functions in the present EMD model, which indicates that this simple model is not 
versatile enough to simultaneously cover in a quantitative way both the magnetic and 
electric sectors of the QGP. 

With these limitations in sight, we expect that the present bottom-up holographic model 
will be mostly useful to understand the effects of magnetic fields on the QGP within the 
range T ~ 150 — 400 MeV and eB < 1 GeV 2 . 

III. HOLOGRAPHIC QCD THERMODYNAMICS AT NONZERO MAGNETIC 
FIELD 

In this Section the results for the holographic equation of state at nonzero magnetic field 
are presented. The formulas needed to compute the observables shown below were presented 
in the last Section. Here, we define the pressure as the temperature integral of the entropy 
density performed while keeping the magnetic field fixed 16 



(25) 


where we took a low reference temperature, T re f = 22 MeV, in agreement with what was 
done in [61, 63] to obtain the fit for the dilaton potential and the gravitational constant (19). 
By doing so, the holographic curves for the pressure in Fig. 4 (and also Fig. 2) actually 
correspond to differences with respect to reference pressures calculated at ^ref for each value 
of the magnetic field. 

In Fig. 4 we show our holographic results for the normalized entropy density, s/T 3 , 
and pressure, p, and compare them to recent lattice data [53] for eB = 0, 0.3, and 0.6 
GeV 2 . It is important to remark, however, that the above convention to calculate the 
16 As discussed in detail in Section 2 of [53] this corresponds to the isotropic pressure in the so-called 
“R-scheme” where the magnetic field is kept fixed during compression. Also, this corresponds to the 
anisotropic pressure in the direction of the magnetic field in the so-called “<E>-scheme” where the magnetic 
flux is kept fixed during compression. 
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pressure is not exactly the same used in [53] since in (25) the pressure (difference) vanishes 
at T = T re f = 22 MeV while in the calculation carried out in [53] the pressure goes like 
~ O (( eB ) 4 ) for T —y 0 and, therefore, one should expect that the differences between these 
two calculations 17 become more pronounced at low temperatures and large magnetic fields, 
as seen in Fig. 4. However, even for eB = 0.6 GeV 2 , we do find a reasonable agreement for 
the pressure at large temperatures (T > 200 MeV). 
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FIG. 4. (Color online) Holographic calculation for the normalized entropy density, s/T' 3 , and 
pressure, p, in the presence of an external magnetic field. The solid, dashed, and dot-dashed curves 
correspond to magnetic fields eB = 0, 0.3, and 0.6 GeV 2 , respectively. The data points correspond 
to the lattice calculations for these quantities performed in [53]. 


On the other hand, when it conies to the ratio s/T 3 , the agreement between our holo¬ 
graphic results and the lattice is only at the qualitative level. This is in part due to the 
uncertainties in the holographic description of this observable already at B = 0: the holo¬ 
graphic model parameters were chosen to describe the lattice data for the pressure and the 
speed of sound squared at B = 0 and not 18 s/T 3 . In any case, one can see that s/T 3 
increases with an increasing magnetic held, which is the general behavior observed on the 
lattice [53]. Moreover, note that the curve s/T 3 becomes steeper near the transition region 

17 Note that in [53] the pressure was obtained from the renormalized free energy density. Here, we could have 
done the analogous holographic procedure by calculating the free energy density from the holographically 
renormalized on-shell action for the EMD model. This is, however, a much more laborious calculation 
than the one we have carried out here where we first calculated the entropy density using the Bekenstein- 

Hawking’s relation (18) and then we calculated the pressure (difference) using Eq. (25). 

18 Probably a better agreement with JS / 0 lattice data may be obtained by improving the choice of the 


model parameters through a global fit to B = CJ^attice data for the pressure, the entropy density, the 
speed of sound, and the trace anomaly. 













for increasing values of the magnetic field, which is again in agreement with the general 
trend observed on the lattice [46]. 


eB [GeV 2 ] 

T c {eB) [MeV] 

0 

158.2 

0.1 

157.6 

0.2 

154.9 

0.3 

153.2 

0.4 

151.3 

0.5 

149.9 


TABLE I. Deconfinement temperature (defined by the inflection point of s/T 3 ) for different values 
of the magnetic field in the bottom-up holographic model. 



FIG. 5. Deconfinement temperature (defined by the inflection point of s/T 3 ) for different values of 
the magnetic field in the bottom-up holographic model. The data points correspond to the lattice 
calculation performed in [53]. 

As discussed in [53], the inflection point of s/T 3 may be used to characterize the crossover 
temperature as a function of the magnetic held 19 . Correspondingly, the peak in T8t(s/T 3 ) 

19 Since the crossover is not a genuine phase transition, the free energy is analytic in the region where the 
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may be used to estimate the crossover temperature as a function of the magnetic field in 
our holographic model. We used our results for s/T 3 to find how the crossover temperature 
changes with a magnetic field and the results are displayed in table I and in Fig. 5. One 
can see in Fig. 5 that in our model the crossover temperature decreases with an increasing 
magnetic field, as found on the lattice [17, 53], but a quantitative agreement with the data 
from [53] occurs only for eB < 0.3 GeV 2 . 

Some general comments regarding the crossover found in our holographic model are in 
order at this point. Depending on the chosen dilaton potential, the black hole solutions 
may or may not have a minimum temperature, as detailed discussed, for instance, in Refs. 
[76, 77, 100, 107]. In the case there is some minimum temperature below which the black 
hole solutions do not exist, the system generally features a first order Hawking-Page phase 
transition [108] to the thermal gas phase at some critical temperature a little bit higher than 
the minimum temperature for the existence of the black hole solutions. Also, in this case, 
the black hole solutions are not unique and there is at least one unstable branch of black hole 
solutions above this minimum temperature. But for some choices of the dilaton potential 
the temperature of the black hole solutions may monotonically decrease as a function of the 
radial position of the horizon until going to zero, in which case the black hole solutions are 
unique and thermodynamically preferred over the thermal gas solution and the system does 
not feature any phase transition at nonzero temperature (at least at zero magnetic field and 
vanishing chemical potentials): this is the case realized in our EMD model. Note also this 
is indeed the adequate situation to mimic the QCD crossover instead of the pure Yang-Mills 
first order phase transition. In fact, by analyzing our dilaton potential according to the 
general criteria discussed in [107], one notes that in the deep infrared our dilaton potential 
goes like V(<f> —> oo) ~ _ e o. 606 <^ j n case a t each finite value of temperature (at zero 

magnetic field and vanishing chemical potentials) there exists a unique black hole solution 
and this corresponds to the true ground state of the system, having a larger pressure than 
the thermal gas solution. Moreover, since within the region of the (T, R)-phase diagram 
analyzed in our manuscript the pressure of the plasma increases with B (as also seen on the 
lattice, see Fig. 4), within this region the black hole solutions are always thermodynamically 
degrees of freedom change from a hadron gas to a deconfined plasma. Thus, the definition of the crossover 
temperature T c depends on the observable one uses to characterize it. Different observables can give in 
principle different values for T c and one may use them to obtain a band defining the crossover region 
[17, 53], 
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preferred and do correspond to the true ground state of the system. 

As a technical detail, in order to obtain the curves in Fig. 4 we used a large grid of 
initial conditions with 720,000 points taking 900 equally spaced points in the (^-direction 
starting from </> 0 = 0.3 and going up to 0 O = 7.8, and 800 equally spaced points in the 
T 5 —direction starting from ~b —— 0 and going up to 73 —— 0.99. A large number 
of points was required to obtain sufficiently smooth curves for s/T 3 that allowed for the 
extraction of the crossover temperature and its dependence on the magnetic held. However, 
smooth curves for p could be obtained using much smaller (and faster) grids. 


IV. CONCLUDING REMARKS AND PERSPECTIVES 

In this paper we developed, for the first time, a bottom-up holographic model that pro¬ 
vides a quantitative description of the crossover behavior observed in the equation of state 
and in the magnetic susceptibility of a QCD plasma with (2 + l)-havors at zero magnetic 
held. We employed this model to study how an Abelian magnetic held B affects the ther¬ 
modynamic properties of this strongly coupled plasma (at zero chemical potentials). In the 
presence of the magnetic held the plasma becomes anisotropic and we used the inflection 
point of the holographically calculated s/T 3 curve to determine how the crossover temper¬ 
ature is affected by the external magnetic held. We found that the crossover temperature 
decreases with an increasing magnetic held, which agrees with the general behavior recently 
observed on the lattice. Our model calculations display some level of quantitative agreement 
with the lattice data for values of the magnetic held up to eB < 0.3 GeV 2 , which is the 
expected range achieved in ultrarelativistic heavy ion collisions. 

We believe that this agreement with the lattice data can be further improved toward 
larger values of eB if one tries to carefully match the lattice thermodynamic calculations at 
B = 0 by simultaneously taking into account different observables such as the pressure and 
the speed of sound squared, as we have done in the present approach, with the addition of 
the entropy density and the trace anomaly in a global ht; in this sense, our choice for the 
holographic model parameters (fixed at B = 0) may be systematically improved. 

An interesting feature of our holographic model that distinguishes it from other con¬ 
structions (such as [109, 110]) is that the suppression of the crossover temperature with 
the external magnetic held found here is directly tied to a quantitative description of near 
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crossover lattice QCD thermodynamics at B = 0. It would be desirable to generalize the 
present holographic model by taking into account the contribution of the chiral condensate. 
Moreover, motivated by the recent studies in Refs. [46, 111], one could also investigate if 
this model indicates the existence of a critical point in the (T, R)-plane at higher values of 
the magnetic held 20 . 

The holographic setup constructed here may be employed to obtain estimates for the mag¬ 
netic held dependence of many other physical observables relevant to the strongly coupled 
QGP. For instance, one could generalize the calculation of transport coefficients performed 
in [71] and obtain a quantitative estimate of how the anisotropic shear (and bulk) viscosity 
coefficients vary with the external magnetic held around the QCD crossover transition. 

Recently, the effects of an external magnetic held on the equilibration dynamics of strongly 
coupled plasmas have been studied using holography [73, 112]. In this context, it would 
be interesting to see how the quasinormal mode spectrum in our nonconformal plasma 
varies with an external magnetic held. Given that our model can capture the nonconformal 
behavior of the QGP near the crossover transition, with and without the external magnetic 
held, a detailed study of the quasinormal modes in this model may shed some light on the 
thermalization process that takes place in an anisotropic nonconformal strongly magnetized 
QGP. We hope to report results in this direction in the near future. 
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Appendix A: Electric susceptibility and conductivity for different coupling functions 


In order to check further limitations of the present EMD model (some of which have 
been discussed at the end of Section IID), we compare in this Appendix the results for 
the magnetic susceptibility, and also the electric susceptibility and DC electric conductivity 
for two different profiles of the Maxwell-Dilaton electric coupling function /(0). The first 
profile is given in Eq. (24), which was fixed by fitting lattice data [52] for the magnetic 
susceptibility at B = 0, as discussed before. The second profile was fixed in Ref. [65] by 
fitting lattice data [113] for the electric susceptibility also at B = 0, 

/(0) = 0.0193 sech(—100 0) + 0.0722 sech(l(T 7 0). (Al) 


At B = 0, the holographic formulas for the electric susceptibility and the DC electric 
conductivity are given respectively by, 

1 s 1 


Q 

X2 


T 2 16v r 2 T 3 /(0) /“ dr e -2ab)/-i(0( r )) ’ 

a Q 27r v / h(p/(0o) 


(A2) 


(A3) 


T K 2 

and we refer the reader to consult Ref. [65] for a discussion on the derivation of these 
formulas 21 . 

One can see from the results shown in Fig 6 that a simple EMD holographic model cannot 
give simultaneously good quantitative descriptions of electric and magnetic field response 
functions: by adjusting the electric coupling /(0) in order to fit the magnetic susceptibility 
at B = 0, one is able to attain a good description of the QCD thermodynamics at finite 
B , as shown in Section III, but response functions to an applied electric field are not well 
described in a quantitative way within such prescription. On the other hand, if one adjusts 
the electric coupling /(0) in order to match the electric susceptibility, one is not able to 
obtain a good quantitative agreement with lattice data for the magnetic susceptibility. It 
would be certainly interesting to think about the construction of some holographic model 
versatile enough to quantitatively cover the entire electric-magnetic sector of the QGP, which 
is something that our simple EMD model is not able to do. We must remark, however, that 
up to now, our EMD model is the only holographic approach available in the literature 


We use that a(rn) = ao = 0 and 0(rjj) = 0o- 
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FIG. 6. (Color online) EMD magnetic susceptibility (top left), electric susceptibility (top right) and 
DC electric conductivity (bottom) for two different choices of the Maxwell-Dilaton electric coupling 
function f(4>): the full curves were obtained by using f(<p) given in Eq. (24), while the dashed 
curves were obtained by employing given in Eq. (Al). All the lattice data displayed in these 
plots refer to (2 + l)-flavor QCD (lattice data for the electric conductivity are taken from [114]). 

which is able to match in a quantitative way the behavior of many magnetic field related 
observables calculated on the lattice. 
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